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Government's  proprietary  interests  in  and  relating  thereto  shall  not  be  Im¬ 
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ABSTRACT 

H 

I 

A  "semi -physical  lit"  to  the  equation  of  state  of  air  &a  generated  from 
<  tabulated  values.  The  iorm  Is 

P  =  (r  -  1)  P  e 
y  -  r  (p,  *) 

p  =  mass/volume 
e  =  energy/maas 

A  listing  of  the  FORTRAN  subroutine  AIR  Is  given,  as  well  as  a  graph 
of  tabulated  and  subroutine-computed  values  of  (y  -  1). 


t 


*» 
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1.  A  SUBROUTINE  FOR  THE  EQUATION  OF  STATE  OF  AIR. 

Computer  calculation*  of  expanding  fireballs  in  the  atmosphere  require 
some  form  of  equation  of  state  for  air  which  is  valid  over  many  orders  of 
magnitude  variation  in  density  and  energy.  Marty  SRELL  calculations  have 
been  made  using  an  equation  of  state  which  appears  much  toe  simple.  How¬ 
ever,  these  calculations  haw  given  results  for  fireball  growth  and  rise 
which  agree  amazingly  well  with  experiments.  Unfortunately,  these  results 
cannot  be  expected  to  give  reliable  pictures  of  temperature  within  the  fire¬ 
ball  simply  because  of  the  inadequacies  in  the  equation  of  state.  The 
distribution  of  temperatures  within  the  fireball  Is  extremely  important  in 
fireball  calculations  since  it  determines  electron  density.  This  note 
describes  an  interim  effort  to  provide  a  good  equation  of  state  for  computer 
codes. 

Previous  SHELL  calculations  have  used  r,  if  e  <  t.  ,  or  r, 

1  ion  *  1 Z 

if  e  >  eioa>  where  is  an  input  value.  The  two  values  of  y  were 

chosen  before  each  problem  by  referring  to  the  real  equations  of  state.  We 
decided  to  try  a  simple  fit  to  (y-1)  for  an  ideal  gas  equation  of  state,  as 
computed  by  MBS.  1  The  variables  which  define  this  equation  of  state  are 
chosen  as  material  density  and  energy  density  (per  mas*)  rather 
density  and  temperature.  Hydrodynamics  calculations  ttlways  involve 
material  and  energy  densities,  and  temperature  is  not  used  as  an  inde¬ 
pendent  variable.  *  ,  . 


The  first  look  at  possible  forrnsfor  tills  tit  was  directed  to  thir  density 
dependence.  Crossplot*  of  (y-l)  vsrsraa  density  for  dilferaml  energies  re¬ 
vealed  an  approximate  behavior  of  the  form  *'*  -  '  :‘vr 


♦This  fact  rules  out  the  use  of  analytic  tits  to  the  equation  of  state  of  air 
using  temperature  aaan  independent  variable  us 

Harold  Brode.  ■  ju*3*T  «■.>  •'  * 
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<r-»  -  [t<p„  - 1|  (£) 


p  =  sea  level  density 

=  1.293  xlO"3  gm/cm3 


The  exponent  a  determines  density  spreading,  and  is  aero  below  e  =  1 
jerk  per  megagram  or  10*®  ergs  per  gram.  Above  e  =  1,  a  je)  is 
roughly  linear  in  log  ( e )  (see  figure  1).  Assuming  the  density  dependence 
discovered,  the  problem  reduces  to  finding  a  fit  for  y  (p^,  e).  The  energy 
dependence  of  fy-1)  is  influenced  by  the  dissociation  of  the  major  con¬ 
stituent  molecules  (only  oxygen  and  nitrogen  contribute  much  to  the 
thermodynamic  properties).  There  are  four  different  "regimes"  which 
are  defined  below: 

I.  Nothing  dissociated  0  <  e  <  *•  5 
II.  Only  oxygen  dissociated  ~  9  <  e  <  *  30 
HI.  Everything  dissociated  ~  50  <  c  <  ~  150 
IV.  Electrons  contributing  ~  150  <  e  <  «• 

Each  "regime**  has  Its  own  energy  dependence,  and  the  thermodynamic 
properties  of  air  are  those  of  whichever  '*  regime"  represents  the  composition 
of  the  air.  Transitions  in. properties  are  made  as  the  composition  changes. 

Regimes  X  and  U  are  affected  primarily  by  the  increasing  occupation  of 
vibrational  and  electronic  levels  of.  mdiaspcictsd  molecules.  The  energy 
dependence  of  regime  HI.  is  due  to  increasing  occupation  of  clpctroadc  energy 
levels,  and  regime  IV  Is  the  threshold  of  the  F ermi-Thoznao  description,  ? 

Thermodynamically,  (y-U  is  the  same  as  2/n*  where  u  is  the  n amber 
of  degrees  of  freedom.  In  a  gas  whose  molecules  or  atoms  have  internal 

iht  13  -.rt-  .sit-y*  ut- •  Vj  i3;  J:  c-  ?•  iici  .*  ,v-»  • 

-~~7  ~  -i*  y.’-.’-^a-scV  tt'  a.*  u-ns.'i^--'  ‘.-.1  i'pxiAv 

*See,  fcr. instance j  "LAolU#*  oh  TFD  thepry *  6iKa;i 

'  ’  -  #  "  vj''.  .  . 
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i  structure,  n  seed  sot  be  an  integer.  For  example,  an  oxygen  or  nitrogen 
molecule  ie  rotating  but  not  vibrating  at  room  temperature,  and  has  an  n  of 
3+2=5.  An  energy  increase  will  make  n  change  from  5  to  7  (although 
i  the  molecule  may  dissociate  before  n  reaches  7).  Regimes  1  and  II  repre¬ 
sent  this  change,  plus  electronic  excitation,  for  0_  and  N,  . 

4U  b 

The  hump  of  increasing  r  at  e  =  8  is  due  to  the  decrease  in  n  for 
oxygen  from  7  to  3  as  the  molecules  dissociate  and  can  no  longer  rotate  or 
vibrate.  A  similar  Increase  in  y  occurs  as  the  nitrogen  dissociates,  and 
again  as  the  electrons*  (with  no  internal  energy  levels)  become  numerous. 

A  fit  So  y  (po>  e)  for  the  entire  interval  could  be  generated  by  a  large 
polynomial,  but  there  is  more  physical  significance  In  a  series  of  experi¬ 
mental  fits  for  each  regime  which  are  "turned  on  and  off*  as  the  composition 
changes. 

In  practice,  "on  and  off  switches"  of  the  form 


1  * 


have  been  found  to  cause  suitable  transitions  at  the  "humps  M  between 
functions  which  have  been  chosen  to  fit  the  data  closely  between  "humps." 
These  switches  have  values  from  0  to  1  with  a  certain-width  (~6Ae) 


centered  about  the  energy  eQ  and  have  no  particular  physical  significance 
in  this  work.  It  seems  reasonable  to  sjsppose,  however,  that  e-  < is'- re¬ 
lated  in  some  way  to  the  dissociation  energy  and  Ae  to  the  dissociation 
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II  W2  choose  to  lit  the  dota  with  exponentials  of  the  form 


Ae 


■tj  c. 


then  oar  method  leads  ns  to  the  following  form  for  the  equation  of  state 


P  «  *  1)  pe 


f  -«/«,  “«/* 2  -*/«,, 

T  -  1  =|a  +  be  1f1+ce  ^  (1-^)  +  de  f2  + 


«{e) 


o  fe)  »  ax  1^  +  a2  (1  -  fj)  {1  -  f^)  +  {t^i 

{We  might  expect  the  density  spread  to  he  different  for  each  regime.) 

-  ,  » 
s  P!  log  *S  «2  =  P2  log  <*3  **  P3. 

The  "switches"  are  given  by 


(Oxygen) 


.1 


\***j&*f+l 

.  £  -T 


(Nltrofep)  ^ 


%?  *X**fC^T 

'  z?>  ■  ««  » 


t 


>  -  J  . 
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{Electrons) 


+  l 


i  This  is  still  not  complicated  enough  to  give  t  close  description,  since  each 
transition  energy  and  width  decreases  with  decreasing  density.  This  density 
dependence  was  represented  by 


e02-»20  ~ 


etc.  and 


Ae02-»*20  ' 


etc. 


The  constants  were  evaluated  by  a  very  straightforward  technique: 
guess  values,  plot  corves,  guess  new  valves.  Toward  die  end,  the  curves 
i  were  competed  by  machine,  but  die  same  technique  was  followed.  After 
many  tedious  attempts,  each  coming  a  little  closer,  it  was  decided  to  stop, 
call  this  an  interim  equation  of  state,  and  continue  with  more  dignified 
techniques.  The  resells  of  the  program  to  date  indicate  that  the  approach 
1  is  worthwhile.  Rattier  than  list  values  for  each  constant  {which  will  prqb- 
■  ably  change  later)  we  list  a  subroutine  (appendix  2}  in  FORTRAN  which  will 
!  produce  fr~D  with  an  error  generally  from  0-2  percent  with  occasional 
increases  to  5  percent.  The  tabulated  and  fitted  values  of  jy-1)  are  given 
in  figure  2.  The  equations  with  numbers  appear  in  appendix  II.  Tempera- 

•v  . 

i  tores  above  1.5  ev  were  not  given  much  cmsiderafica,  since  radiation  Crcns- 


port  Influences  things  above  that  temperature,  and  our  planted  application 
is  for  T’simple"  hydrodynamics  only.  It  may  be  possible  to  continue  the 
description  of  regime  Vf  to  Inflate  energy^  wpgjg aTFD  oqastioa  ol  sdnt 

Farmer  cpplicafions  of  fids  fit  axe  B^tm^Ms.^'W^t  this  equation  of 
state,  -one. can  use  SB^Lfst  fireball  ccjciputi^S^  sad'  prdutiy  believe 
The  same  form  h^rheiE^^/oAer.gases  and  trapes. 
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Lt.  J.  D.  Garcia  ha*  laggested  that  altering  the  energies  of  transition  may 
account  for  nonequllibrium  effects  in  dissociation  as  the  fireball  expands 
rapidly.  This  seems  quite  possible,  once  it  can  be  determined  how  to 
compute  the  new  transition  energy.  Techniques  to  determine  the  coefficients 
from  theoretical  methods  (as  opposed  to  curve-fitting)  may  discover  new 
freedom  in  generating  equations  of  state  for  other  gases. 

Finally,  the  continuation  of  this  work  to  include  equations  of  state  for 
earth,  fire,  and  water  will  allow  calculations  for  all  possible  substances. 
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A  LASTING  OF  THE  FORTRAN  SUBROUTINE  AIR 
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Appendix  II 

The  equation  with  the  present  numbers: 

T-1  =  {.161  *  .255  e-''4- “  «!  •  Zaoc-'^U-J,** 

.lS7.-e/M’Slz*  .050  f3} 
a(e)  =  ,>84S  loglQ  *  +  .032(3 


e(e) 


c  -  e 


e  Aet  «  i 


£_2_  C 


*  Ae2  +  1 


y  a  - 

f2>  I°810  *  +  .1 

Ei  = 

8.5  + 

,357  1o510 

„  .*0S 

AEj  = 

.975 

(3 

.0157 

e2  = 

45. 

(t) 

„  .085 

Ae2  = 

4.0 

Ik) 

*2 


fill 


e  -Ac 


3 


e3  =  160. 
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